Analysis of nucleus-nucleus collisions at high energies and Random Matrix Theory 
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We propose a novel statistical approach to the analysis of experimental data obtained in nucleus- 
nucleus collisions at high energies which borrows from methods developed within the context of 
Random Matrix Theory. It is applied to the detection of correlations in a system of secondary 
particles. We find good agreement between the results obtained in this way and a standard analysis 
based on the method of effective mass spectra and two-pair correlation function often used in high 
energy physics. The method introduced here is free from unwanted background contributions. 
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I. INTRODUCTION 

There is currently an enormous effort underway to 
detect signals of possible transitions between different 
phases of strongly interacting matter produced in high 
energy nucleus- nucleus collisions (of [lljQ). It is expected 
that in central collisions, at energies that are and will 
be soon available at CERN Super Proton Synchrotron 
(SPS), BNL Relativistic Heavy Ion CoUider (RHIC), and 
CERN Large Hadron Collider (LHC), the nuclear density 
may exceed by a few orders of magnitude the density of 
stable nuclei. At such extreme conditions one would ex- 
pect that a final product of heavy ion collisions could 
present a composite system that consists of free hadrons, 
quarks and gluons. However, quarks and gluons cannot 
be unambiguously identified with present detectors. In 
addition, the identification, for example, of a quark-gluon 
plasma is hindered due to a multiplicity of secondary par- 
ticles created at these collisions. In fact, there are numer- 
ous additional mechanisms of particle creation that mask 
the presence of a quark-gluon plasma. A natural ques- 
tion arises as to how to identify a useful signal that would 
be unambiguously associated with a certain physical pro- 
cess? 

Amongst the most popular traditional methods of 
analysing data produced at high-energy nucleus- nucleus 
collisions are: i) the correlation analysis Q , ii) the analy- 
sis of missing masses [4j| and effective mass spectra Q , iii) 
the interference method of identical particles p etc. It is 
currently believed that by measuring event-by-event fluc- 
tuations it would be possible to observe anomalies from 
the onset of deconfinement [3]. Often though, results 
based on such methods are sensitive to assumptions made 
concerning the background measurements and mecha- 
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nisms included in the corresponding model. For exam- 
ple, one of the popular approaches to studying particle 
production in heavy ion collisions is through statistical 
mechanics techniques (cf (alsli^)- However, recent the- 
oretical analyses of the multiplicity fluctuations demon- 
strate that corresponding results are different in different 
ensembles and it is sensitive to conservation laws obeyed 
by a statistical system |ll| . 

Transport theory has played an important role in the 
interpretation of experimental results in heavy ion re- 
actions at various energy domains (see, for example, 
RefsUallSl). The underlying concept of transport models 
is based on numerical solution of the Boltzmann equation 
(the Vlasov-Uehling-Uhlenbeck, the Boltzmann-Uehling- 
Uhlenbeck (BUU) etc). However, exact numerical solu- 
tions, for example, of the BUU equation are very difficult 
to obtain and, therefore, various approximate treatments 
have been developed [Ij, [Q- Up to now there are no 
stringent proofs that such approximations really yield a 
solution of the full transport problem. The only proof 
is that numerical solutions provide the asymptotic limit 
at equilibrium, which can be deduced from the Boltz- 
mann equation (the BUU equation without potential and 
Pauli-blocking factor) (see details in Refll2l). We recall 
that the BUU equation is a generalization of the Fokker- 
Planck equation associated with the Brownian motion of 
particles in a classical Coulomb gas model. Dyson dis- 
covered that this model can be described well in terms of 
the Gaussian ensembles of the Random Matrix Theory 
(RMT) fl4| . This fact implies that one might expect that 
the RMT could be useful for analysis of experimental re- 
sults obtained in heavy ion reactions at high energies. 

The RMT was originally introduced to explain the sta- 
tistical fluctuations of neutron resonances in compound 
nuclei [l^l (see also Ref. [l6( ). There the precise heavy- 
compound-nuclear Hamiltonian is unknown or rather 
poorly known, and there is a large number of strongly 
interacting degrees of freedom. Wigner first suggested re- 
placing it by an ensemble of Hamiltonians which describe 



all possible interactions [17|. The theory assumes that 
the Hamiltonian belongs to an ensemble of random matri- 
ces that are consistent with the fundamental symmetries 
of the system. In particular, since the nuclear interaction 
preserves time-reversal symmetry, the relevant ensemble 
is the Gaussian Orthogonal Ensemble (GOE). Whereas if 
time-reversal symmetry were broken, the Gaussian Uni- 
tary Ensemble (GUE) would be the relevant ensemble. 
The GOE and GUE correspond to ensembles of real sym- 
metric matrices and of complex Hermitian matrices, re- 
spectively. 

It is noteworthy that during the last twenty years the 
RMT grew into the powerful new statistical theory of 
fluctuations in a variety of physical problems. For ex- 
ample, the RMT has become a very successful theoret- 
ical method for analysing fluctuation properties found 
in the data from atoms and nuclei, quantum dots, and 
many other systems (see, for example, Ref 1 181) . In various 
fields, the Dyson-Mehta statistical measures are the most 
often used to quantify a system's correlations and to de- 
termine what information the fluctuations contain. They 
can be used with any discrete data to search for correla- 
tions. These measures do not depend on the background 
of measurements and used in the context of RMT give 
universal forms depending only on the fundamental sym- 
metries preserved. Their only requirement is that local 
mean densities (or secular behaviors) be understood and 
their effects be removed. Furthermore, a change of fluctu- 
ation properties of a system under consideration, induced 
by a change its symmetry properties, can be detected by 
the RMT tools unambiguously (see below). Therefore, 
these tools could provide a way of detecting the transi- 
tion between different phases (with different fluctuation 
properties) of strongly interacting matter produced in 
high energy nucleus- nucleus collisions. 

The main aim of the present paper is to introduce 
the Dyson-Mehta statistical measures for analysis of ex- 
perimental data from nucleus- nucleus collisions at high 
energies in order to detect effects produced by corre- 
lations. More specifically, to demonstrate that they 
are very sensitive to spacing correlations present in the 
nucleus-nucleus collision data, and serve a useful purpose 
of independent complementary approach to the standard 
tools designed for such an analysis in high energy physics. 
Note, in general, the standard tools are based on different 
concepts, in contrast to the Dyson-Mehta measures that 
are formulated within unique, mathematically rigorous 
theory. In Section II we review how to detect the mani- 
festation of correlations with the aid of nearest neighbor 
spacing distribution (NND) from the data obtained in 
light nuclei collisions in Dubna experiments. To check the 
validity of our findings by dint of the NND in Section III 
we analyse the same data within the method of effective 
mass widely used in the analysis of data from nucleus- 
nucleus collisions at high energies. Section IV is devoted 
to analysis of density-density correlations based on the 
Dyson-Mehta statistical measures and the two-pair cor- 
relation function often used in high energy physics. The 



main conclusions are summarized in Section V. The pre- 
liminary results of our analysis were presented in [19| . 



II. NEAREST-NEIGHBOR SPACING 
DISTRIBUTION 

A. Some experimental details 

To demonstrate the feasibility of using the Dyson- 
Mehta statistical measures in the context of high en- 
ergy nuclear collisions, we make use of the experimental 
data that have been obtained with the 2-in propane bub- 
ble chamber of High Energy Physics Laboratory (LHE) , 
JINR [20l . \2u . [23 |. The chamber, placed in a magnetic 
field of 1.5 T, was exposed to beams of light relativistic 
carbon nuclei at the Dubna Synchrophasotron. Nearly 
all secondary particles, emitted in a An total solid angle, 
were detected in the chamber. 
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FIG. 1: The momentuni distribution of the secondary charge 
particles. 

All negatively charged particles, except those identified 
as electrons, are considered as 7r~-mesons. The contami- 
nation from misidentified electrons and negative strange 
particles does not exceed 5% and 1%, respectively. On 
FiglT] the momentum distribution dN/d\p\ of the sec- 
ondary charged particles is displayed. The maximum of 
the particle production is observed at \p\ ~ 0.5 — 0.7 
GeV/c. The particle momenta were calculated from the 
particle trajectories in a magnetic field taking into ac- 
count ionization and radiation losses. In particular, the 
uncertainty in the momentum value were estimated tak- 
ing account of the effects of multiple Coulomb scattering 
and bremsstrahlung radiation. The average uncertainty 
in the momentum and the angle measurements varies 
as: {Ap/p) = (11.5 ± 0.3)%, {A9) - 0.8". The mini- 
mum momentum for pion registration is about 70 MeV/c 
in lab frame (we consider all kinematic variables in the 
lab frame). The protons were selected by a statistical 
method applied to all positive particles with a momen- 



turn of \p\ > 500 MeV/c (slow protons with \p\ < 700 
MeV/c were identified by ionization in the chamber). 

In this experiment, there are 37792 -'^^C+^'^^C interac- 
tion events at a momentum of 4.2A GeV/c (for greater 
discussion of the details see [2ll, [22, [231) containing 
7740 events with more than ten tracks of charged par- 
ticles. The separation method of ^^C-|-"'^^C collisions 
in propane {C^H^), data processing, identification of 
particles and discussion of corrections are described in 
detail in [H. About 70% of i2c-|-i2c events have 
been identified with the aid of this method. The set 
residue of events was separated statistically between car- 
bon-|-carbon and carbon-|-hydrogen collisions with the 
aid of the event "weight". The "weights" were estimated 
from known cross sections for inelastic carbon+carbon 
and carbon-|- hydrogen collisions. 



B. Basic remarks 

If one supposes that the momenta of secondary parti- 
cles produced in nucleus-nucleus collisions may be treated 
in analogy with eigenstates of a composite system, just 
like the eigenstates of the compound nucleus, then the 
statistical analysis methods ,15] introduced by Dyson and 
Mehta can be applied to the LHE collision data. By anal- 
ogy with compound nuclei discussed in the Introduction, 
one may assume that the observed spectra (Fig(T]) be- 
long to the effective QCD Hamiltonian which describes 
the composite system, but is unknown. We point out that 
the secondary particle momenta are well determined in 
the experiment. It is, therefore, natural to use the mo- 
mentum as a proper variable for our analysis in order 
to avoid possible errors at the transformation from the 
momentum to the energy, which requires an accurate de- 
termination of the corresponding mass value. 

Based on this supposition, the ordered sequence of "en- 
ergy levels" {Ei},i = 1,...,N comes from the ordered 
sequence of the detected momenta of the secondary par- 
ticles. Hereafter, for the sake of convenience we associate 
a symbol " E" with an absolute value of the momentum 
\p\. To analyse the statistical fiuctuations of the spec- 
trum {Ei}, it is necessary to separate its smoothed aver- 
age part whose behaviour is nonuniversal and cannot be 
described by the RMT 25]. A separation of fluctuations 
of a quantum spectrum can be based on the analysis of 
the density of states below some threshold E 
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p{E) =Y.^{E - E,) . 



(1) 



We can define a staircase function 



S{E)^ f p{E')dE' = Y,e{E - E'), (2) 

giving the number of points on the momentum axis which 



are below or equal to E{pr \p\). Here 

e{x) = 



for X <Q 
for a; > 1 



(3) 



We separate S{E) in a smooth part (^{E) and the re- 
minder that will define the fiuctuating part Sfi{E) 



SiE) = ({E) + SaiE), 



(4) 



such that the integral of S{i{E) is zero. Here, the smooth 
function C{E) gives the mean number of eigenvalues less 
than or equal to E of the exact eigenvalue distribution 
S{E). The smooth part ({E) can be determined ei- 
ther from semiclassical arguments or by using a polyno- 
mial/spline interpolation for the exact staircase function. 
From this sequence a new one is obtained by the unfold- 
ing procedure of the original spectrum {Ei} through the 
mapping E ^ x (see details in Sec lII Cp 



aE^), 



t = l. 



,N. 



(5) 



The effect of the mapping is that the sequence {xi} has 
on average a constant mean spacing (or a constant den- 
sity), irrespective of the particular form of the function 
C{E) [25|. Note that there is no any combinatorics in- 
volved in a such procedure. That which remains in the 
sequence are the fiuctuations away from unit mean. Next, 
one defines the spacing Si = x^+i — Xi between two adja- 
cent points. Collecting Si in a histogram, one obtains the 
probability density or the NND. We stress that, in gen- 
eral, this procedure does not involve any uncertainty or 
spurious contributions and deals with a direct processing 
of physical data. 

Bohigas et al [26[ conjectured that the RMT describes 
the statistical fluctuations of a quantum system whose 
classical dynamics is chaotic. Quantum spectra of such 
systems manifest a strong repulsion (anticrossing) be- 
tween quantum levels, while in non-chaotic (regular) sys- 
tems crossings are a dominant feature of spectra. In 
turn, the crossings are usually observed where there is no 
mixing between states that are characterized by different 
good quantum numbers, and the anticrossings signal a 
strong mixing due to a perturbation brought about by 
either external or internal sources (cf Refs[l5l[l6l). 

Thus, if the "events" {xi} are independent, i.e., cor- 
relations in the system under consideration are absent, 
the form of the histogram must follow p{s) = exp(— s) 
known as the Poisson density. On the other hand, if the 
levels are repelled (anticorrelated) as in the GOE, the 
density is approximately given by the Wigner surmise 
form p{s) — ^sexp(— js^). In other words, any cor- 
relations that produces the deviation from the regular 
pattern (Poisson distribution): production of a collec- 
tive state (resonance), or some structural changes in the 
system under consideration would be uniquely identified 
from the change of the histogram shape. In particular, 
the transition from one probability density to the other 
has been used in nuclear structure physics to study the 
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FIG. 2: Nearest-neighbor spacing momentum distribution (histogram) p(s) for different regions of measured momenta: (a) 
and (b) correspond to 0.15 < \p\ < 1.14 GeV/c; (c) and (d) correspond to 1.14 < \p\ < 4.0 GeV/c; (e) and (f) correspond to 
4.0 < IpI < 7.5 GeV/c. The total distribution (it includes data related to all charged secondary particles) consists in (a), (c), 
and (e) panels. The distribution for proton momenta consists in (b), (d), and (f) panels. The Poisson and the Wigner surmise 
distributions are connected by dot-dashed and solid lines, respectively. 



stabilization of octupolc deformed shapes and transition 
from the chaotic to the regular pattern in the classical 
limit ,27j . Therefore, such an analysis can provide the 
first hint of some structural changes at different param- 
eters of the system under consideration, in particular, in 
different energy (momentum) intervals. 



C. The experimental NND histograms 

In order to reveal structural changes in the momentum 
distribution we determine the smooth part C(^) using a 
polynomial function of the fifth order to interpolate the 
exact staircase function (see also Eq.((3])): 



c(i^) = ^afei?^ 



(6) 



k=0 



where E = \p\ are the experimental momentum values 
in the given event. We recall that momenta are ordered 
in ascending series. The parameters ak were optimized 



with the aid of the program MINUIT. Next, we obtain 
a new spectra {xi} by the unfolding procedure of the 
original spectrum {Ei} through the mapping _B — > a; by 
means of the polynomial function ^ with the optimized 
coefficients 0^: Xi — C{Ei),i = 1, ...,N. Having a set of 
variables {xi}, one is able to determine the set of spac- 
ings {si}. Since all events are independent we apply the 
same procedure for the other events and obtain the new 
independent sets of spacings {si}. To reduce the uncer- 
tainty in the data we put an additional constraint: the 
distributions of various spacings Si from the 7740 events 
satisfy the condition of x^ per degree of freedom is less 
than unity (in our case six parameters). Collecting all 
independent spacing together we require the fulfillment 
of the following conditions: 



I p{s)ds = 1, / sp{s)ds = 1. 
^0 



(7) 



The first condition provides the normalization of the 
probability, while the second one gives the normalization 
of the mean distance to unity. 



To clearly recognize correlations we divided the total 
set of spacings {si} on three sets, in correspondence with 
three regions of the measured momenta: a)0.1 < \p\ < 
1.14 GeV/c (region I); b)1.14 < \p\ < 4.0 GeV/c (re- 
gion II); c)4.0 < \p\ < 7.5 GeV/c (region III) (see FigH]). 
The region boundaries were determined with the require- 
ment that the shape of the spacing density p{s) does not 
change in the region under consideration. For example, 
the decrease of the upper bound of the region I does not 
change the Poisson distribution (see below). However, 
the increase of this boundary leads to the deviation from 
the latter distribution. The same criteria has been ap- 
plied for the other two regions. Note that there is no a 
prescribed procedure how to define such regions. How- 
ever, the empirical approach described above proved to 
be useful in data processing for various systems at the 
RMT analysis (see references in RefsUa. UM . We also 
separated the distributions associated with the proton 
momenta in order to illuminate possible effects related 
to the mass difference between proton and pions. Al- 
though there is the distinction between the total and 
proton NNDs, a general pattern is similar in the both 
cases (compare left and right columns of Figl2). Thus, 
the mass difference does not change the overall conclu- 
sions and it is important for the utility of the approach. 

The NNDs in the region 0.15-1.14 GeV/c are displayed 
on Figs. la, b, where the minimum value of the proton 
momentum is 0.15 GeV/c. One observes an uncorrelated 
pattern: the function p{s) has the Poisson distribution. 
In this region the momentum distribution was defined 
with a high accuracy. The region II covers the values 
1.14-4.0 GeV/c (Figs.lc,d). This region corresponds to 
the intermediate situation, where the spacing distribu- 
tion lies between the Poisson and the Wigner distribu- 
tions. One may conclude that in this region there is al- 
ready the onset of some correlations. These correlations 
manifest themselves much stronger in the region 4.0-7.5 
GeV/c (region III, see Figs.le,f). The region HI is char- 
acterized by the Wigner-type distribution for the spacing 
probability. 

From the above analysis one may conclude that in this 
particular experiment the method enables us to reveal 
structural changes in the system at particular energy re- 
gions. At high energies (the third interval) the particles 
are produced at relatively small spatial separation, since 
Ap is large and one can make use of the uncertainty re- 
lation ApAx ~ h. Thus, one might expect enhanced 
correlations. These correlations continue to persist at 
the second interval but with a lesser intensity. And, fi- 
nally, in the first interval the secondary particles move 
independently in the spatial separation defined by the 
kinematics. To illuminate the physical nature of this ob- 
servation we consider below the method of effective mass, 
traditionally used to process the data produced at high 
energy nucleus-nucleus collisions. 

To complete the discussion we calculate the integrated 
momentum spacing density for experimental data in or- 
der to confirm that all distributions have on average a 




FIG. 3: (Color online) Cumulative spacing momentum distri- 
butions for different regions of measured momenta: a) 0.1 < 
IpI < 1.14 GeV/c (solid line); b)1.14 < \p\ < 4.0GeV/c 
(dashed line); c)4.0 < \p\ < 7.5 GeV/c (dot-dashed hne). 
The solid line nearly coincides with the result for the Poisson 
distribution, while the dot-dashed line is close to the Wigner 
surmise distribution. 



constant mean spacing (=1). This is an additional test 
which verifies the convergence of our calculations. Fig[2] 
shows the integrated density, demonstrating a gradual 
transition from a Poisson-like density toward a Wigner- 
like density with the increase of the absolute value of the 
momentum distributions. 



III. METHOD OF EFFECTIVE MASS 

In this section we consider the method of effective mass 
(MEM) [5| , which is a standard tool to extract the infor- 
mation on the correlation between secondary particles. 

Before we proceed a few remarks are in order: in vari- 
ous theoretical approaches it is usually assumed that at 
high-energy nucleus-nucleus collisions: i) a majority of 
produced secondary pions are emitted basically through 
the mechanism of production and decay of the light res- 
onances; ii) a significant portion of the protons are pro- 
duced as a result of A-isobar decays |28i|. Further, it is 
also assumed that such processes consist of two steps. 
First, (n — k) particles and the resonance are produced. 
Second, the resonance decays on fc-particles and one may 
expect that due to kinematics there are some correla- 
tions between these /c-particles. In order to extract these 
correlations we are forced to consider all possible combi- 
nations of the fc-particle-participants and compare with 
the known resonance masses. 

In particular, in the MEM one considers the effective 
mass for k secondary particles which create a cluster (a 
resonance) as 



fc ~. 1/2 



^ 1=1 1=1 ■' 



(8) 



If the condition Afj, = MJl'^'' holds, one may conclude that 



a resonance with the mass M™** is identified in the data. 
In addition, one assumes that each identified resonance 
contributes to the total cross section with its own weight 
(probabiUty) . To carry on this idea, each resonance is 
approximated by the Breit-Wigner distribution with the 
identified mass and the resonance width. Varying the 
weights of identified resonances, one is aiming to reach 
the best agreement with the observed inclusive cross sec- 
tions of resonances production. This procedure proved 
to be reasonable at small multiplicity of the secondary 
particles or at relatively small energies per particle (< 10 
GeV/c). The situation changes drastically at the large 
multiplicity, where number of the secondary particles 
is more than ten in each event. For example, suppose 
that we have N secondary particles in each event. The 
most favorable case is where each pair of the secondary 
particles creates a resonance; the number of resonances 
Nres = N/2 in the event. However, we have to consider 
all independent combinations between pairs in the event 
and, therefore, the total number of physical and "spuri- 
ous" resonances is Ntot = N*{N — l)/2. It results in the 
ratio r = Nres /Ntot ~ 0.1 (a "useful" signal) for ten par- 
ticles in the event. Evidently, this procedure gives rise to 
a large spurious contribution to the analysis. 
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FIG. 4: (Color online) The distributions (crosses) of (yr^p)-, 
{tt~p)-, (pp)-, and (7r"'"7r~)-pairs emitted in the interval of 
the momentum distribution of the secondary particles 0.1 < 
\p\ < 1.14 GeV/c. The background is shown by solid line. 
The total number of pairs is normalized to unity. 



pairs constructed from the particles produced in differ- 
ent events only. Note, that this procedure may lead to 
the violation of the energy and momentum law conser- 
vation due to the mixing of different events. Finally, the 
resonance production is associated with the existence of 
the difference between the distribution of all independent 
combinations {Ntot = N * {N — l)/2) in the event and 
the background distribution for the pairs. 

Keeping this in mind, we consider the distributions of 
charged particle pair s emitted in ^^C-|-^^C -interactions 
at 4.2 A GeV/c [H ^. Figs HIS] demonstrate the dis- 
tributions of (tt+p)-, {Tr~p)-, ipp)-, and (7r+7r~)-pairs 
emitted in three ranges of the momentum distribution 
of the secondary particles: 0.1-1.14 (FigH]); 1.14 - 4.0 
(FigE]) and 4.0-7.5 GeV/c (FigE]). AU distributions are 
normalized to the total number of pairs. Note, the low- 
est threshold for the effective masses scale is determined 
by the individual mass of each participant and their mo- 
menta, according to Eq.(I5]). 

One observes that in the interval of momentum 0.1 — 
1.14 GeV/c (FigU) no clear-cut distinction exists be- 
tween experimental and background distributions. In 
this interval we have obtained very good statistical con- 
ditions: 41615 (tt+p)-, 43626 (tt-^)-, 52992 {pp)- and 
39112 (7r+7r~)-pairs. One concludes that there is no a 
manifestation of a resonance production. We recall that 
in this interval the RMT approach produces the Pois- 
son distribution for the behavior of the p(s)-function (see 
Figl2^). In other words, the NND indicates on the ab- 
sence of any correlations in this energy interval. 
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Further, to extract the information on the resonance 
production from data it is necessary to evaluate a back- 
ground contribution. It is indeed a difficult task which 
cannot be solved completely. Therefore, to construct the 
background the method of mixing events [29| is usually 
used. In this method the background is defined by the 



FIG. 5; (Color online) Similar to Fig|4]for the interval 1.14 — 
4.0 GeV/c. 

In the second interval (see FiglS]) we have obtained 173 
(tt+p)-, 16470 (TT-p)-, 167094 {pp)- and only 16 (tt+tt-)- 
pairs (are not shown). There is some deviation of the 



experimental distribution from the background only for 
(7r^p)-pairs. Although the error bars are noticeable, 
there is a clear indication on the resonance production 
at 1.1 < Meff < 2 (GeV/c^). We believe that it is 
connected with the production of the A+^-isobars with 
masses 771^++ = 1.232 and 1.650 (GeV/c^). Note that the 
observation of A++-isobars with masses toa++ = 1.232 
GeV/c^ has been reported in [2^, while the ones with 
masses 1.650 GeV/c^ are observed for the first time in 
the present data. The statistics are not enough, how- 
ever, to support strongly the latter observation. There 
are no visible deviations between the experimental data 
and the background for the distribution of the {tt~p)- 
and (pp)-pairs. In this interval the RMT approach pro- 
duces a visible deviation from the Poisson distribution for 
the behavior of p(s)-function (see FigO:;), i.e., the NND 
signals upon the onset of correlations. The agreement 
between predictions obtained in two different approaches 
confirms that the NND is able to provide a hint of the 
resonance production. 

In the third interval of the momentum distribution of 
charged secondary particles we have 10 {■K~p)- (are not 
shown) and 9522 (pp)-pairs (FigE]). Here, the (tt+p)- 
pairs are absent. One observes a clear cut distinction 
between the background and experimental distribution of 
(pp)-pairs in 4.0 — 7.5 GeV/c. The deviation of the signal 
relative to the background is above 20%. However, there 
is no a solid basis to associate such a str ong deviation 
with a production of di-baryon resonances [30j , since the 
inclusive cross section of such "resonances" would exceed 
essentially those that are predicted by various theoretical 
models. 
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FIG. 6: (Color online) Similar to Fig|4]for (pp)-pairs emitted 
in the interval of the momentum distribution of the secondary 
particles 4.0 - 7.5 GeV/c. 

It is well known that in this interval the stripping pro- 
tons are the dominant ones (with a small contribution of 
deuterons, tritons and others) |22| . These protons carry 
a maximum momentum near the value 4.2 GeV/c. It 
results in very small deviations of the particle trajecto- 



ries in the magnetic field of the setup. In fact, it is the 
worst situation for the accurate determination of the er- 
rors in the momentum distribution. The RMT approach 
produces in this interval a distribution of the density p(s) 
close to the Wigner surmise form (see Figl2j;). As stressed 
above, such a distribution is associated with the break- 
ing of regularity in the spectral properties of a quantum 
system due to either external or internal sources. In [19| 
we have already mentioned that the onset of the Wigner 
distribution for the density p{s) (breaking of the regu- 
larity) could indicate the presence of errors in the mea- 
surements, i.e., the correlations introduced by external 
perturbations. 
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FIG. 7: Left: The proton distributions dN/d\p\ as a function 
of the measured momentum in three different regions: 0.1 — 
1.14 (a); 1.14-4.0 (c) and 4.0-7.5 (e) (GeV/c). Right: The 
proton distributions dN/dO as a function of the angle in three 
different regions: 0.1 - 1.14 (b); 1.14 - 4.0 (d) and 4.0 - 7.5 
(f) (GeV/c). The angular distribution is normalized to unity. 

To gain a better insight into the nature of correlations, 
it is worthy of note that in all considered distributions 
there is an evident dominance of the pp-pairs. To trace 
the evolution of the pp correlations we select only the mo- 
mentum and angular distributions of the protons in three 
intervals (see Figl?]) . In the first interval the angular dis- 
tribution (see FiglTb) of the pairs covers almost all angles 
of the semisphere with some concentration around ^ 50". 
We recall that the NND (see Figll}]) displays the Poisson 
distribution in this region, which is typical for a nonin- 
teracting system. The absence of the correlations for 
lower-energy proton pairs is due to Fermi protons (with 
the average momentum ~ 300 MeV/c) which move as 
independent particles. The angular distribution is deter- 
mined by the kinematics (the momentum and the energy 
conservation laws). In the second interval the momen- 
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turn distribution of the pairs (Fig[7l;) is spread smoothly 
over all considered values of the momentum (1.14 — 4.0 
GeV/c) . There is some concentration of the emitted pairs 
in the angular distribution, which covers a solid angle 
~ 20° (Fig[7ti). The protons can be divided on partic- 
ipants and spectators. The participants are due to the 
A-isobar and excited nucleon decay, and also consist of 
the protons which acquire momenta in direct nuclear re- 
actions. Although there are some signature of spectators, 
they are dominant in the third region. For the second re- 
gion we obtain the NND pattern which is neither the 
Poisson nor the Wigncr type (FigOi). In the third inter- 
val (Figsl7]3,f) one observes that stripping protons have 
similar momenta and almost zero angle in the distribu- 
tion. Evidently, under such conditions, one may expect 
a large probability for the interaction in the final state, 
which leads to the narrow peak appearances in the effec- 
tive mass spectrum of the proton pairs. Such interaction 
effects in a final state are well known for the particle 
production and decay process at high energies [3l|. We 
recall that for this interval the RMT approach provides 
the Wigner distribution (Fig[2f) for the behavior of the 
density p{s). In other words, the NND indicates that in 
this region there are strong correlations that break the 
regularity in the spectra. The maximum of the corre- 
lation is observed to increase with the increasing of the 
kinetic energy. This may indicate that the energetic pro- 
tons create clusters due to strong pairing effects at short 
time intervals. We will return to this point in the next 
section. 

The comparison of the NNDs with the MEM analy- 
sis manifests in fact that there are evident correlations 
between behavior of the density p(s) in different energy 
(momentum) intervals and the appearance of new sources 
that break the regularity in the momentum distribution 
of the charged particles. The NNDs discussed in the pre- 
vious section indicate the onset of correlations (a pro- 
duction of resonances) at the intermediate region unam- 
biguously, via the deviation of the spacing distribution 
from the Poisson limit (see FiglJ). According to the 
MEM analysis the stripping protons are responsible for 
the strong correlations in the third interval. We stress 
that the presence of strong correlations is predicted in- 
dependently by the NND which exhibits the Wigner dis- 
tribution in this region. Presently, however, we cannot 
distinguish between correlations introduced by the phys- 
ical process and by possible errors in the measurements. 
This is one of the main challenges for the RMT approach 
in high energy physics. 



IV. CORRELATIONS 



In order to characterize the degree of correlations for 
a stationary spectrum with unit average spacing Dyson 



introduced the k-level correlation functions 
Rk{xi,...,Xk) = (Af^fc)! / ■■■! Pk{xi, ...,XN)dXk+l-.dxN 

1 <k<N, (9) 

where Pis!{xi, ...,xiq)dxi...dxj^ gives the probability of 
having one eigenvalue at xi, another at X2■■■^ another at 
xjq each within the interval {xi, Xi+dxi\. By integrating 
Pfc(xi, ..., xn) over all variables but one, in the limit N — > 
oo, one obtains the ensemble averaged density 



P{xi) = 



Pk{xi, ...,XN)dx2...dx 



N 



(10) 



which is normalized to unity. From Eq.Q it follows 
that i?i(a;i) = Np{xi) and Rk{xi, ..., Xk)dxi....dxk is the 
probability, irrespective of the index, of finding one level 
within of each of the intervals [xi , Xi + ds] . From the 
above definition it follows that Ri{x) — 1. With the aid 
of the definition ([9]), by integrating Rk+i one obtains 



Rk+iixi,...,Xk+i)dxk+i = [N - k)Rk{xi,...,Xk) 



(11) 




FIG. 8: (Color online) Number variance E (L) for three dif- 
ferent regions of the experimental spacing momentum distri- 
butions: a)0.1 < IpI < 1.14 GeV/c (solid line); b)1.14 < \p\ < 
4.0GeV/c (dashed line); c)4.0 < \p\ < 7.5 GeV/c (dot-dashed 
line) . 

It is difficult to work directly with the Rk functions. 
One important and more convenient measure of correla- 
tion that was introduced is based on the number statistic 
n{E, L) which is defined to be the number of levels in an 
energy interval of length L and centered at the energy 
(momentum) E: 

r.B+L/2 

p{E)dE. (12) 

-L/2 

Since the spectrum was unfolded (see Eq.®), the ensem- 
ble average number statistic {n{E, L)) = L is indepen- 
dent of the spectrum. However, the variance of n(_E, L) 

Y?{L)^{[n{E,L)-{n{E,L))Y), (13) 



n{E,L)= f 

J E 



does depend on the spectrum considered. For the Poisson 
density (see [ij]) 

Y?{L)^L, (14) 

and for the GOE one the exact asymptotic expression is 



s^W = 6, 



ln{2TTL) + 7 + 1 + -[Si{nL)]'^ - -SiinL) 



' cos(27rL) - Ci{2TTL) + tt'^L 



1 - -Si{2TTL) 



(15) 



Here 7 is the Euler constant and Si, Ci are the sine 
and cosine integrals, respectively. The number variance 
S^(L) calculated using the optimal implementation of 
the definition in Eq. p^ is shown in Fig. [51 Indeed, this 
quantity manifests the Poisson statistics (Eq. fM)) ) for ex- 
perimental spectra with a low momentum distribution. 
On the other hand, one again observes a clear indication 
of the presence of correlations for large momenta. 



A 




FIG. 9: (Color online) Two-point correlation function R2{s = 
L) for different regions of measured momenta: a)0.1 < |p| < 
1.14 GeV/c (solid line); b)1.14 < \p\ < 4.0GeV/c (dashed 
line); c)4.0 < \p\ < 7.5 GeV/c (dot-dashed line). The solid 
lines, denoted as P and GOE, display the characteristic limits 
for Poisson and GOE ensembles, respectively. 

For the analysis of fluctuations it is more convenient 
to use pure k-point functions [32| 



RkiL) 



Rkixi,...,Xk)dxi...dxk 



(16) 



The function Rk{L)/kl gives the probability that an in- 
terval of length L (for small L) contains k levels. In the 
RMT most emphasis has been put on the two-point corre- 
lation function i?2 (2^1,2:2) or density-density correlation 
function. The two-point correlation function is the prob- 
ability density to find two eigenvalues Xi and Xj at two 
given energies (momenta), irrespective of the position of 
all other eigenvalues. The function R2{xi,X2) depends 
only on the relative variables s = xi — X2 defined in 



Secini The variance S^ is connected to R2 through the 
following relation 



R2{L) - 



/ R2{xi,X2)dxidx2 ^ 2 / (L — s)i?2(s)ds 
"'0 Jo 



Y.^L) + L{L - 1) . 



The two-level correlation function R2{xi,X2) determines 
the basic fluctuation measures related to Wigner's level 
repulsion and the Dyson-Mehta long-range order, i.e., 
larg e correlations between distant levels. Bohigas et al 
[331 provided a thorough analysis of level repulsion and 
long-range correlations (rigidity) for different correlation 
functions. To understand the distinct role played by level 
repulsion and long-range order in the momentum density, 
we compare our numerical results with analytical expres- 
sions from Table 1 of Ref 32 for the Poisson ensembles 
(there is neither level repulsion nor long-range order) and 
for the GOE case (this ensemble exhibits both features). 
The two-point correlation function R2{xi,X2) calculated 
from Eq. pT)) is shown in Fig[5] Even though there are 
small deviations from the Poisson {R2(L) = L^) and the 
GOE {R2{L) = 7t'^L^/18) predictions, the experimental 
results for the momentum distributions reproduce sur- 
prisingly well both limits. 

We recall that the analysis is based on the unfolded 
spectra obtained from the experimental momentum dis- 
tribution (Sec. II). In contrast to the MEM, the above 
procedure does not introduce any " spurious" correlations 
due to combinatorics. In addition, each event is treated 
separately. However, the two-point correlation function 
R2{xi,X2) clearly exhibits the presence of correlations in 
the experimental data. It is another beautiful example 
how the correlations can show up in the Dyson-Mehta 
statistics. Thus, the transition from the Poisson limit to 
the GOE implies the onset of strong correlations at the 
region HI. 

As was stressed above the two-point correlation func- 
tion i?2(s = L) indicates upon large correlations between 
distant levels (momentum). It appears that the tran- 
sition to the GOE limit might signal on the formation 
of clusters. These clusters should consist of energetic 
protons mainly, which seems move pairwise in the third 
interval (see Fig[7f). Note, however, that the NNDs for 
the third region indicate on the Wigner repulsion. The 
physical nature of this phenomenon may be understood 
with the aid of the Hanbury-Brown-Twiss (HBT) anal- 
ysis of the proton correlations [33]. The HBT analysis 
focuses on the spatial distribution of the radiation area, 
being very sensitive, however, to the statistical assump- 
tions upon the background distribution. Evidently, such 
analysis represents a separate problem itself and is be- 
yond the scope of the present paper. 

Nevertheless, it is interesting to note, that the standard 
pair-correlation function (see, for example, RefiSJI and 
references therein) 



R{Vi,V2) = cr 



d^(j/dyidy2 
{da / dyi){d(7 / dy2) 



1, 



(18) 
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used for the analysis of data in high-energy physics, might 
be useful to compare with the RMT results. Here, the 
quantity a is the cross section of the inclusive reaction 
and y — ^In ^_p^ is the rapidity, which depends on the 
particle energy E and its longitudinal momentum P| | . 
The rapidity is one of the main characteristics widely 
used in relativistic nuclear physics (see 0, Q ) . In partic- 
ular, the change of the reference frame leads to a trivial 
shift Ay in the rapidity. The pair-correlation function 
manifests the difference between the probability density 
of two-particle events and the product of the probabil- 
ity densities of independent particle events. It vanishes 
if the particle rapidities are independent, i.e., the corre- 
lations are absent. In a way, this function is similar in 
spirit to the two-point correlation function R2{xi, X2), or 
density-density correlation function of the RMT. There- 
fore, it is useful to compare the predictions obtained by 
dint of the RMT tools and with the aid of the standard 
pair-correlation function. 

R 
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yi 

FIG. 10: Integrated two-pair correlation functions for parti- 
cles obtained in ^^C-|-^^C- interactions (see the text) for differ- 
ent regions of measured momenta: a)0.1 < \p\ < 1.14 GeV/c; 
b)1.14 < IpI < 4.0GeV/c; c)4.0 < \p\ < 7.5 GeV/c. 

Fig [To] displays results for the standard pair-correlation 
function, obtained for ^^C-|-^^C-interactions. The func- 
tion R{yi,y2) was integrated over one of the variables, 
say yi, and we consider the dependence of this func- 
tion on 2/2- For different momentum distributions, there 
are three intervals of integration for the variable j/i: 
a)for 0.1 < IpI < 1.14 GeV/c the function R{yi,y2) 
is integrated in the interval —1.0 < yi < 2.5; b)for 
1.14 < IpI < 4.0 GeV/c it is integrated in the interval 
0.1 < yi < 2.4; and c)for 4.0 < |p| < 7.5 GeV/c it is 
integrated in the interval 2.5 < yi < 3.5. The results 
for the function R = J Riyi,y2)dyi clearly indicate the 
presence of correlations between particles in the region 
4.0 < IpI < 7.5 GeV/c : there is a strong deviation from 



zero in the interval 2.7 < 2/2 < 3.2. While this function 
does not exhibit any signature of correlation in the second 
region, one observes a strong correlations between rapidi- 
ties in the third region. We recall that in this interval the 
NND exhibits the Wigner distribution (short-range cor- 
rellations), typical for strongly correlated system. Thus, 
we speculate that in the region III there is a formation 
of the short-living proton pairs (see also SeclIIip. 
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V. SUMMARY 

Using the Dyson-Mehta statistical measures for the 
first time for analysis of data in high energy physics, we 
have studied the experimental data from ^^C-|-^^C col- 
lisions at 4.2 A GeV/c. The high accuracy of the data 
provides a reliable basis for the analysis of the momen- 
tum distribution of the secondary charged particles pro- 
duced in this reaction. The NND analysis (Sec. II) in- 
dicates on the onset of correlations in the momentum 
distribution interval 1.14 < |p| < 4.0 GeV/c (region II). 
The results of calculations for different correlations func- 
tions (Sees. II, IV) evidently demonstrate the presence of 
strong correlations in the momentum distribution inter- 
val 4.0 < IpI < 7.5 GeV/c (region III). 

We have analysed the same data with the aid of the 
method of effective mass, one of the standard tools for 
analysis of data produced at nucleus-nucleus collisions 
at high energies (see Sec lIIip . The results of this anal- 
ysis exhibits the presence of strong correlations, inde- 
pendently found by dint of the RMT approach. In the 
region II we found the production of well known A"*""*"- 
isobars with masses m/^++ = 1.232 and 1.650 (GeV/c^). 
The dominance of the proton pairs with zero angle in the 
momentum distribution interval 4.0 < |p| < 7.5 GeV/c 
could be attributed to the interaction effects in the final 
state. It appears that these effects may lead to the for- 
mation of the short-living proton pairs. Note that the 
method of effective mass involves, however, some addi- 
tional assumptions that induce unavoidable uncertainty 
in the analysis due to a large nonphysical contribution 
to the data. This may cause additional errors, which are 
difhcult to exclude with the aid of the approaches widely 
used for the analysis of nucleus- nucleus collisions at high 
energies. In contrast, the RMT approach is independent 
of such spurious contributions. The only requirement of 
the RMT is that the local mean densities (or secular be- 
haviors) should be understood and their effects removed 
through the use of the unfolding procedure (see Sec ITT)) . 
which is a routine procedure for the RMT. By dint of rel- 
atively simple and transparent mathematical tools, the 
RMT approach detects the appearance of new sources 
that break the regularity in the momentum distribution. 

Furthermore, the predictions based on the RMT analy- 
sis are consistent with the predictions based on the stan- 
dard pair-correlation function R (see Sec lIVp . which is 
another tool for the analysis of data produced at nucleus- 
nucleus collisions at high energies. In fact, the RMT 
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two-point correlation function magnifies the presence of 
correlations manifested in the standard pair-correlation 
function (compare FigslQl andflOh). We stress that the 
uncertainty of the results encountered in the theoreti- 
cal analyses of the multiplicity fluctuations in different 
statistical ensembles [ll| does not appear in the RMT 
approach at all. The Dyson-Mehta statistical measures 
rely heavily on the fundamental symmetries preserved in 
the nucleus- nucleus collisions. 

We conclude that the RMT approach is free from var- 
ious assumptions concerning the background of the mea- 
surements and it provides reliable information about cor- 
relations induced by external or internal perturbations. 
We recall that it is expected that at central heavy-ion 
collisions at relativistic energies the multiplicity of sec- 
ondary particles will increase drastically. Evidently, the 
standard tools like the MEM approach (Sec lIIip would 
encounter a large uncertainty due to the increase of spu- 
rious background contributions in the analysis of data. 
In contrast, the efficiency of the RMT approach increases 
with the increase of the multiplicity of the secondary par- 
ticles in each event taken individually. We stress that 



the RMT approach is based on a unique, mathematically 
rigorous theory at the limit of A'^ — > oo. Therefore, we 
believe that our analysis provides the basis for further de- 
velopment of quantitative methods based on the ideas of 
the RMT for data processing produced at nucleus- nucleus 
collisions at high energies. 
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